Riemann normal coordinates, Fermi reference system and the geodesic deviation 

equation 



Alexander I. Nesterov * 

We obtain the integral formulae for computing the tetrads and metric components in Riemann 
normal coordinates and Fermi coordinate system of an observer in arbitrary motion. Our approach 
admits essential enlarging the range of validity of these coordinates. The results obtained are applied 
to the geodesic deviation in the field of a weak plane gravitational wave and the computation of 
plane- wave metric in Fermi normal coordinates. 
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I. INTRODUCTION 



The well known Riemann normal coordinates satisfy the conditions g^lp = an( i = at a point p, the 

origin of the coordinate system. In these cordinates the metric is presented in powers of the canonical parameter 
of the geodesies out from p. In 1922 Fermi 0] showed that it is possible to generalize the Riemann coordinates in 
such a way that the Christoffcll symbols vanish along given any curve in a Riemannian manifold, leaving the metric 
there rectangular. For a geodesic curve Misner and Manasse introduced a Fermi normal coordinates, which satisfy 
Fermi conditions along a given geodesic and calculated the hrst-order expansion of the connection coefficients and 
the second-order expansion of the metric in powers of proper distance normal to the geodesic. Their construction 
is a special case of Fermi coordinates defined by Synge ||. The last ones form a natural coordinate system for a 
nonrotating accelerated observer. 

In 1973 Misner, Thorn and Wheeler (MTW) § defined a Fermi coordinates for an accelerated and rotating observer 
and calculated the first-order expansion of the metric. Extending MTW work Mitskievich and Nesterov || and later Ni 
and Zimmerman [[| obtained the second-order expansion of the metric and the first-order expansions of the connection 
coefficients. Li and Ni derived the third-order expansion of the metric and the second-order expansions of the 
connection coefficients. 

The size of neighbourhood V, where Riemann and Fermi normal coordinates do not involve singularaties, is deter- 
mined by conditions 



u <C min 



o 

RajS/3 

~o 



RcfjSfl I 1 / 2 | Ra^Sfl,^ 



(I 



where u is a canonical parameter along the geodesies from the origin of coordinates. The first condition u <C \ Ra-yg/3 
\~ x / 2 determines the size of V where the curvature has not yet caused geodesies to cross each other. The second 
condition determines the domain where the curvature does not change essentially. For instance, for gravitational 
waves with wavelength A the Riemann tensor is ~ A e~xjp(ikx) / \ 2 , where A is the dimensionless amplitude. This 
yields 

u < minjA/VA, A}. 

Generally it is assumed A < 10 -18 . This means that the size of V is restricted by uo <C A. So the application of 
Riemann (or Fermi) coordinates to the modern experiments may be very restrictive since A is often supposed being in 
the order of 300 km. Thus for enlarging the range of validity by a factor 1/yA (which is about 10 9 in our example) 
it is necesary to take into account derivatives of any order of the Riemann tensor. Some years ago Marzlin || has 
considered Fermi coordinates for weak gravitational field <? p „ = r/^ + h^, {h^l <C 1 and dervived the metric as a 
Taylor expansion valid to all orders in the geodesic distance from the world line of Fermi observer. 

Here we develop a new approach to Riemann and Fermi coordinates, which admits enlarging the range of validity 

of these coordinates up to u <C | R ai s/3 | -1 ^ 2 . The point is to use the integral formulae. Recently this approach has 
been used to calculate a geometric phase shift for a light beam propagating in the field of a weak gravitational wave 

The paper is organized as follows. In Section 2 we obtain the integral formulae for computing the tetrad and metric 
components in Riemann normal and Fermi coordinates. In Section 3 we compute the plane-wave metric in Fermi 
normal coordinates and analyse the geodesic deviation equation for this case. 

In this paper we use the space-time signature (+, — , — , — ); Greek indices run from to 3, Latin a, 6, c from 1 to 2 
and k from 1 to 3. 



II. GENERAL RESULTS 



A. Riemann normal coordinates 



For arbitrary point po in some neighbourhood V(po) there is the unique geodesic 7(u) connecting po and p which, 
using exponential mapping, we may write as p = exp 7 ( u ) (it£), where £ = d/du — £ Q e( Q ), u being the canonical 
parameter while basis is parallely propagated along "f(u). The Riemann normal coordinates with the origin at the 
point po, are defined as X a = w£ Q , where u£ = exp~^(p). The following basic equations summarize the most 
important properties of the Riemannian normal coordinates: 
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3(a) (Po) = (d/dX a ) P0 = {Ajd/dxt*),,, A J = (dxP/dX a ) P0 (1) 

r a /37 = o, r a 7 x^ = o, r^ 7 = -| i Q (/37)5 + • • • , (2) 

lo s 1 4 

= + - RafS/3 X 1 X + — R al S/3,p X 1 X X^ 

I/ 16 ° \ r ^ 

T . 6 i? Q75/3 a, M +^ R ai s P R p x^ X~<X S X X X^ + ■■■ (3) 



5! 

The size of neighbourhood V, where these coordinates do not involve singularaties, is determined by conditions 

^ ■ J 1 I I 1 (A s 

"«™ 1 — L± — > ■ (4) 

[ I R ai S0 1 1/2 \R ai Sf3,n\) 

The first condition u <C | Ra-ySp | -1 / 2 determines the size of V where the curvature has not yet caused geodesies to 
cross each other. The second condition determines the domain where the curvature does not change essentially. 

For studing Riemann and Fermi normal coordinates one useful technique is the geodesic deviation equation pj. 
Further we apply it for our computation. Let us consider a point po and a frame {e( a )(po)} in this point. For 
arbitrary geodesic 7(11) out from po the equation of geodesic deviation has the following form 

V 2 r/+R(^)£ = 0, (5) 

where £ = d/du is tangent vector to 7(14), rj is vector of deviation and R is operator of curvature. We suppose that 
the basis e( a ) is parallely propagated along *f(u). In tetrad components the equation (Q) takes the form 



= R (a \m,)(^ w ^ h) v {5) (6) 



du 2 

where 7/"-* = rj^e^, etc. 

In Riemann normal coordinates there are natural deviation vectors r\ a = uS^d/dX^ satisfying Eq.(||) p|. Then (j|) 
yields 

J2 (a) „ 1 («) 

We solve this equation by succesive approximations |fl2f . Let us write (0) in the equivalent integral form 



e 



(a) 



= S$ + - dr dT'r'e^R\ S0 C^, (8) 
u Jo Jo 



where we assume el = Sg at the origin of Riemann coordinates. Now we start with 



To improve this first approximation we feed back into the integral (^) getting 



„(") 

e l/3 



= 51 + - T dr f T dT'T'R a lS8 ^^, 
' uJo Jo 



(9) 



Repeteang this process by substituting the new efg back into Eq. (||), we find our solution el""* as infinite series 



E (»)(«) , lrVl 
e/3 , (10) 



in i _ \^ (™^°) 



[ U dr f drV y{ a) &wC?, 
Jo Jo 

+i)(«) 1 r r , , pA f7t5 

e/3 =— dr drr e x R 75/3 £ . 
u Jo Jo 



(2) (a) 1 

e/3 = - 
u 

(n+l)(a) 1 
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(n)(a) ... 

Noting that the order of magnitude of functions ep is given in the neighbourhood V by 

(»)(«). \R alS ^X S \ n mSLX 

60 ~ (2n+l)!! ' 
we conclude that the radius of convergence of series ([l0|) is determined by 

f 1 1 



U 



1^7-5/3 I max 



Now using (f7T0) one obtains for the components of orthonormal tetrad in Riemann coordinates the following 
integral formula: 



- f U dr f dT>T'R a lS ^£, s + 0(R 2 ) 
u Jo Jo 



= 5 + -l dT I dT' T 'R a lS ^e + 0{R 2 ). (11) 



Remark 2.1. Integrating by parts one can write (llll) as 



e {a) - g 



PU PT c\ PU pT pT 

dr dr'R a 7 s P C^-- dr dr' dr" R a lSf) ? + 0(R 2 ). (12) 
Jo Jo u Jo Jo Jo 



This form is more convenient for the comparison with the results obtained in Fermi coordinates (see Subsection B). 
For the covariant derivative of the tetrad we obtain the following useful formula (see Appendix A) 



V ax e M " = -- [ R^xp^rdT + 0(R 2 ). (13) 
u Jo 



With the aid of (|ll|), it is straightforward to derive the metric in Riemann coordinates by using g a p = r)^ v ea 



The result is 



g a p = + - [ dr f dr'r'R a ^C^ + 0(R 2 ), (14) 
u Jo Jo 



and can be rewritten by using Taylor expansion as 

1 ~ s 1 

^ Ra~/8/3 XX + — Rofyop,ft -- ■"■ i gj 



lo s l o s 6o s x 

9a/3 = Va/3 + r Ra~f5f3 X 1 X + — R a7 S0,n X 1 X X^ + — RaySp^,^ X 1 X X X^ H 



that agrees with Eq.(3). 

In fact, we are interested in expressions which should involve (in the right-hand sides) the arbitrary chosen initial 
non-Riemann coordinates only. Then 

i \ 1 f u f T f)X a r)T a Bt^ Bt 5 

4 - s "> + d dT l "'#~&mw&mefr + ° <*■). (», 

The transformation from arbitrary coordinate system to Riemann one takes the form x^ 1 = x fl (po) + A^X V + 0(F) 
(or x^ = x 1 * (po) + A^fu + 0{T) ) and the equations @, (0) can be written as 



el Q) =8$ + -[ dr I dr'r' R\ Sa A x la A JA/ ''ApTf + o(R), (16) 
u Jo Jo 

g a p = Va p + - [ U dr f dr' r' R XlSa A* A^A/A^Ce + o(R), (17) 
u Jo Jo 

where R X jS/3 — R^-ySaix^ipo) + A^ v t'\ while the integration is performed. 
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B. Fermi coordinates 



Let us consider the proper system of reference e( a ) of a single observer which moves along a worldline 7 with an 
acceleration and rotation Hi 



De {a) 



n-e (a)) (18) 



ds 
where 

W = G"t u - G v t" + E^t^ujs 

where E pv ^ s is the axial tensor of Levi-Civita and t p is the unit tangent vector to the line 7; G p being the four- 
acceleration of observer, ui s the four-rotation of the tetrad. Fermi coordinates are defined as X p = 5qS + 8 p £ l u (see 
e.g. Here s is proper time along the observer's world line 7, £ is a unit vector on 7 tangent to the spacelike 

geodesic starting from 7, ortogonal to it and parametrized by proper length u. In fact, £ describes the direction in 
which such a space-like geodesic goes, and it possesses only spatial components different from zero (the temporal 
coordinate is directed along 7). In Fermi coordinates, the covariant components of the corresponding orthonormal 
tetrad parallel transported along the spacelike geodesic, are represented as expansions (see Appendix B) 



e M) = 6% + tt\X l + X - RP ij0 X l X> + i R 0ij0 W k X l X>X k + ± R» ij0 , k X l X>X k + ■■■, 

e M , = 5 p+\ Wito x * xj + R °»p ^ k X l XiX k + -L R» ijp , k X i X 3 X k + • • • , (19) 


where quantities of the type of Q are taken on 7. The expansion up to 3d order of the metric and up to 1st order of 
the connection coefficients (Christoffel symbols) is given by |^-|^] 

a k k a i ' k 

<?oo = 1 + 2QopX p + n ak Q a P X X p + Ro k pQ X X p + R a ijo £l a k X X 3 X 



lOpSL i-SiafeSi pA W* -f JriOkpO -A j\ ' -|- JriaijC. 
1 .,,10 

- Roijo Qo k X z X J X + — 



"g Roijo Slo k X l X 3 X k + — Roijo,k X l X°X 
2 , , 1 „ „*.„,„„, 1 



floi — ^ip^ p + - Rikio X k X l + — Rikio QomX k X l X m + — Rokio Qi m X k X l X m 
3 4 6 

1 u r 1 u 1 

+ 7 -ftifcim " nA A A +T KiklO,m AAA + • • • , 

6 4 
lo , lo , , , , 

ffij — %' + g + ^(ftfciO %m-X" A A'™+ i?jfc;o W A' A'" 1 ) 

+ g Riklj,m X X X" 1 + ■ • • , 

rg„ = ^ + (<^s> 4 + 6° v w x n\ - n v n^+ r p vi0 )x* + ■ • • , 

2 



r « = -o^fe)fe^ +•••• (20) 



For applications it is convenient to write the expansion for the metric tensor in the following form 





.goo = (1 + GXf - (u> X Xf+ R 0kp0 X k X p + R l ij0 e lkp u p X*X 3 X k 
40 . 1 k 

- i?0ij0 G k X l X 3 X + - RoijO,k X t X 3 X 

2o fc lo , , 1 

7 Rokpi X X p + — Ri k io G m X X X + — 
3 4 6 



goi — — (<*> X X)j + - Ro kp i X h X p + — Rikio G m X k X l X m + — Rokia €i mn uj n X k X l X m 

1 u 1 1 



-- R m kli £mnpU P X k X 1 X n + — RiklO.m X k X 1 X* 

6 4 

1 , „„ 1 , _„,.„,„„, 



ffij — %' + - Rikpj X k X P + — ( _Rifcio £jmriW™A' fc A''A' m + i^jfeiO tinin^™ X k X 1 X m ) 

1 , , 

+ g Riklj,m X X X" 1 + • • • , 
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where the expression (w X X)i denotes ey/jO^X and GX := —GiXK 

The radius of convergence of series ( f20|) is determined by (see Introduction) 



ti<iii) = min 



o 

11 1 | RajSp 



|G| ' lEa^l 1 / 2 ' lAryW 



Let us compute the integral formula similar to ([ll]) . In Fermi coordinates the deviation vectors are given by j^j 



The straightforward calculation leads to 



4 Q) = § p + - I" dT f dr'r'e^R x ijp ee, (21) 
u Jo Jo 

o Q) = ^ + + f dr f dr'e^R\ jo ee. (22) 



where 



The iterative scheme (hOh yields 



j ( Q ) 

o(a) Q a e M 



(a) \ ^ (")(«) (°)( Q ) jo , r0n« ci /oo^ 



n=0 



{n tf [a) = 1 f f drY fl^fV, (24) 



'o Jo 

and the radius of convergence of series (p3|) is determined by conditions 



In the first approximation one obtains 

1 f u 



e («) =6? + - dr dr'r'(St + fifl"^^?? + 0(R 2 ), (25) 



e 



u Jo Jo 

(a) 



^=^ + 0°,-^+ / dr / dr'(5S , + <5Sn Q i fO^i*0^e fc + O(iZ 2 ). (26) 



Integrating by parts we combine (|25| ), (|2q ) into one expression 

Jo Jo 



2 



dr dr' I d^'^+Sl^iCr'^^e + OiR 2 ), (27) 
o Jo 



and the componets of metric tensor are given by 



dr / dr'i^x + 5°x%ieT')R x jk 0¥S k 



/ 

cy pU pT pT 

- dr dr' dT' , FJTi fi x + #Slf li ?T' , )& Jkp e i $ k 
u Jo Jo Jo 



2 r 



o Jo 



dr / dT'(^ x +5 a x n^eT')R x jka ¥e 



0{R 2 ), 



dr \ dr' dT"6 p JVaX + 5ln ai eT")R x jkp ?Z k 



This integral formula agrees with the series ( p0| ) . 

The transformation from arbitrary coordinate system {x 11 } to Fermi one takes the form 

a" = ^ l (X°) + Af(X°)X l + ---, 

and in the coordinates {x^} one can write Eq.(|27j) as 

e ( f =S%+ SjWiA^Cu + dr dr'(A- la + AA^-^O^VOi^A/Afc^ 

Jo Jo 

-- ( U dr f dr' f T dT''6 p JA- la +Ax A~ la n\Cr'')R x SKp A/A k K e^ + o(R), 
u Jo Jo Jo 

where it is used the same notation R a /3jg for the curvature tensor in the coordinate system x 11 and 



R' 



iiuX = R\„x (xP(X°) + Af (X°)er") 



while the integration is performed. 



(28) 



III. WEAK PLANE- WAVE METRIC IN FERMI NORMAL COORDINATES 



We compute the plane-wave metric in Fermi normal coordinates for a geodesic observer in the weak gravitational- 
wave field, whose metric tensor is usually written in synchronous coordinates, 

ds 2 = r] lll ydx fl dx l/ + h ao dx a dx b , 

where a and b run from 1 to 2 while 

h a b = h ab (t - z), h 22 = -hn = h + , hi2 = h x . 
Using the definition of the Riemann tensor 

RfiuXcr — ~^(hvX,fi,a ~t~ h[i(7,v lambda h>fiX,v : a ^i/ct,/i,a)i 

one finds that in the linearized theory non-zero components of Riemann tensor are 

RzabZ — RoabO — ~R?,abO = -^h an , 
R^22v = —Rfillu, Rfil2v = Rp,2\v, 

dot being derivative with respect to t. 

The next step is to choose an orthonormal frame along the world line 7 of observer. Taking into account the 
transformation to the Fermi normal coordinates x^ = + 0(h ao ) and that the timelike base vector must be the 
tangent d/dt, one obtains 

e (0) | 7 = d/dT = d/dt, e (1) | 7 = d/dX = d/dx, 
e( 2 ) | 7 = d/dY = d/dy, e (3) | 7 = d/dZ = d/dz, 

where X° = T, X 1 = X, X 2 = Y, X 3 = Z are scttcd. 

Using (E8|) we compute the metric tensor in the Fermi normal coordinates. The result is 
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where 



goo — 1 + H ab X a X b , 

9ab = Vab + F ab Z 2 , 
90a = ^ (^ab + ^ab) X b Z, 

9oz = -\{Hab + F ab )X a X\ 

g 33 = -l + F ab X a X b , 
9a3 — —F ab X b Z 1 



H ab = ^ (h ab (T -Z)- h ab (T) + Zh ab {T)) , 



(29) 
(30) 

(31) 

(32) 

(33) 
(34) 



F a b Hab 



H ab (T,Y)Y 2 dY. 



It is easy to see that H' ab = -F ab , H ab (T,0) = 3F afc (T,0) = (l/2)h ab (T) and 

F a b = 7;({Ftab — F ab )Z)' , 



here and below dot being derivative with respect to T and prime derivative with respect to Z. The nonvanishing 
components of Christoffcl symbols are found to be 



T^a 1 rjf ry2 pa tt vb ^ / tt ttI \ v6 ry 

L 0b-2 H ab A ) L m -tlabX - -{Hab - H ab )X A 

T^a r. Tp vb i Tpl yb ry pia TP V ^ T?' 7^ 

1 33 — ^ab A +-^ab A Z ; 1 36 — ~-^ab Z ~ TJ a 6 

ra rr JTl ry xrfe pO pi3 TT vb 

03 — -H ab JL — £l ab ZjJ>s. , 1 Qa — 1 Qa — tl ab ^ , 

Tab = tF ab Z + -F' ab Z 2 , = -F a bX b , 



1 00 



1 03 



2 

r3 ^ rV \r a x^b 

oo — 2 at A > 

r3 ^ rjV y/) 

03 — —tlnh^ y\ , 



1 33 



r 3 

1 33 



1 p' V a 

2 ab A A ' 



and the computation of the non-zero components of Ricmann tensors in Fermi normal coordinates yields 

-R3a63 = -RoahO = ~ -R3a60 



This leads to 



— 2~{FabZ 2 + 2H ab - (H a b + F a b)Z}. 



RzabZ = RoabO = -RzabQ = -^ lah ^ ~ 



(35) 

(36) 
(37) 
(38) 

(39) 

(40) 

(41) 
(42) 



Thus in Fermi coordinates the curvature is the same function (l/2)ft, a h, but depending on (T — Z). 

Plane morwehromaric gravitational wave. The only nonzero components h a b = A a bsin[k(t — z)\ of wave amplitudes 
are 



Computation yields 



H a b — k A ab 
Fab = jcos(fcT) 



/in = -h-22 = A + sin[/c(i — z)], 
hi2 = h 2 i = A x sin[/c(i - z)]. 



nrT1 ,kZ -sm(kZ) . . JfT1 ,cos(fcZ)-l 
cos(fcT) ——-k + sin(fcT) 



(kzy 



[kZf 



kZ - sin(fcZ) _ cos(fcZ) - 1 + (kZ) 2 /2 



(kZ) 2 
+ sin(fcT) 



(kZf 

cos(fcZ) - 1 fcZ-sin(fcZ) 



(kZf 



(kZf 



(43) 



(44) 
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This leads to 



where we set 



-k 2 A ab X b Z \ cos(fcT) 



.900 - 1 + k 2 A ab X a X b 
.903 = k 2 A ab X a X b (cos(fcT) 



90a 



gab = Vab + k 2 A ab Z 2 \ cos(fcT) 



9a3 



,, m . kZ — s'm(kZ) ., „. cos(fcZ) — 1 

cos(fcT) — — ^ ^+sin(A:T) v y 



(fcZ) 2 v ; (fcZ) 2 

sin(fcZ) - fcZ cos(fcZ) - 1 + (kZ) 2 /2 



(kZ) 2 {kZf 

sin(fcZ) - kZ cos(fcZ) - 1 



sin(fcZ) - kZ cos(fcZ) - 1 + (kZ) 2 /2 



- sin(fcT) 



(fcZ) 2 (kZf 

sin(/cZ) - fcZ cos(fcZ) - 1 



(kzy 



(kZf 



kZ-sin(kZ) 2 cos(kZ) - 1 + (kZ) 2 /2 



(kZ) 2 
sin(kT) 



(kZ) 3 

cos(fcZ) - 1 2 &Z - sin(fcZ) 



(fcZ) S 



-/c 2 A a6 X & Z <^ cos(fcT) 



fcZ-sin(fcZ) 2 cos(fcZ) - 1 + (kZ) 2 /2 



(kZ) 2 
sin(kT) 



(kZf 

cos(kZ) - 1 /cZ-sin(fcZ) 



333 = -1 + k 2 A ab X a X" i cos(fcT) 



fcZ-sin(fcZ) 2 cos(fcZ) - 1 + (kZ) 2 /2 



{kzy 

sin(fcT) 



(kZ) 3 

cos(fcZ) - 1 { 2 kZ- sin(fcZ) 



,4 ab X b = A+(6lX - 5 2 a Y) + A x (SlY + S 2 a X), 
A ab X a X b = A+(X 2 - Y 2 ) + 2A X XY. 



(45) 



(46) 



(47) 



(48) 



(49) 



(50) 



Our results agree with the metric obtained in |13|. 

Geodesic deviation. Let us consider the geodesic deviation of two neighbouring geodesies in the gravitational field 
of a plane weak wave. The separation vector rj satisfies a geodesic deviation equation (||), which reads 



d\ 2 



(51) 



A being a canonical parameter along the geodesic and £ a tangent vector to it. The solution of this equation is given 
by 



[ dr [ dT'R a PlS C P C rfi, +0{R 2 ) 
Jo Jo 



(52) 



where we set t]q = 77° (0) and integration being performed along the basic geodesic line canonically parametrized by A 
The separation magnitude is found to be 



l = lo- 



dr I dT'R aM5 rifao 



(53) 



Now let us consider the geodesic deviation of two neighbouring particles A and B assuming that the the origin 

of Fermi reference frame attached to A's geodesic and for T = the particles are in the plane Z = 0. In this case 



the components of the separation vector are nothing but the coordinates of particle B: r/ a = X a , C /3 = 8q and the 
canonical parameter A being a proper time r along the geodesic line of the observer. Then (153]) reads 
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1 ' A 



l = l °~tJo ^ I dT ' R °-boX a X b . (54) 



(i 



This yields 



l = lo- r h ab(T)X a X b (55) 

'0 



and for the plane monochromatic gravitational wave one obtains the well known result 



I = l - — (A + cos 2tp + A x sin 2tp) sin(fcT). (56) 
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APPENDIX A 

Here we shall obtain the formula (Eq.([lj|) in the text) 

Vd,e {v f = — / rdrR^xs f 5 +o(R). 
u Jo 

Let us start with the Taylor expansion for tetrad in the neighbourhood V(po) 

, o v „ o v 

u being the canonical parameter along geodesies starting in the point po. Using in the Riemann normal coordinates 
the equation of parallel propagation 



de {v f 
du 

we rewrite (1571) as 



Kx^)^ X = 0, (58) 



e {v f =8i-% r^ ff , 7 X°X~< - i T" vtr , y , s X°X">X S + ■■■ + 0(R 2 ). (59) 



The expansion of the connection coefficients is given by 

K x =T" uXii a 1 + i r" tfAi<|I x l x l + ■■■. (60) 

Applying Eq. (jH^), ( |60| ) we obtain the following series for the covariant derivatives of tetrad 
_ o i 

v<9A e o) M =r M „A, CT x° '— r M iy(A (T ) x a + — (r^A,^ - r At l/ (A >CT , 7 ))x°'A" 7 

+^(r"^,«r, 7 ,« - r» u(x ^ s) )x°x^x 5 + o(i? 2 ). (6i) 



Now using the relations 



10 o o 

2 o oo 

- R^ u<j\^x c 'x 1 — (r^^A.CT^ — r M i,(A, C r,7))^ CT x 1 , 

30 o-<5° ° S 

7 -R^ j/crA, 7 <5^ CT ^ 7 ^ = (r M yA : cr,7 : (S — r M 1/ (A, (Ti7! «5))A <T A 7 A , etc 
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one can write the expansion mWj as 



n 

y R^ ^ 7i a, 72 ,..., 7 „^ 71 ^ 72 • • • X^" + ■ ■ ■ + 0(R 2 



(n + 1) 

It is convenient to present this series in the form 

o 



V 8x e (l/) " - V ^ ( , } , , , p- + 0(i? 2 ) (62) 



du n (n + 2)(n!) 

Straightforward calculation yields the following integral representation of Eq. ( |6l| ) 



- f U R^xaCrdr + OiR 2 ). 
u Jo 



Va A e (l/) " - -- / R^ux^rdr + 0{R 2 ). (63) 



Integrating by parts it found to be 

V dx e {l/ f = -[ R^x«CdT + - f dr [ dr' R^xa? + 0(R 2 ). 
Jo u Jo Jo 



(64) 



APPENDIX B 

Here we shall obtain the series (Eq.([l9|) in the text) 

e£° = 6% + fl^X 1 + i R» ij0 X l X> + 1 R 0ij0 W k X l X>X k + 1 fl" ii0|fc + • • • , 

e p " } = <5£ + i fl% p X 4 ^' + i?% P ^ fe X i X^'X fe + ^ R» ijp , k X l X>X k + ■■■. (65) 
We start with the Taylor expansion for tetrad in a world tube surrounding the world line 7 of Fermi observer: 



w being the canonical parameter along spacelike geodesies ortogonal to 7. Using in Fermi coordinates the equation of 
parallel propagation 



^- " I^eWf = 0, (67) 
we rewrite Eq.(|66j) as 

is if °, \ . 

e£° - rv X 1 + - I r% M + r A ^r v A/ 1 x l x l 

I (o v \ 0^ 0^ o^ °x °v °» \ I 

+ 3! \ ^ Vfi i,i,p + r /ii^r^Ap + r fj,iT^xi, P + r M ir CT Air i/ <7P I x l x x p + ■. (68) 

For computing of the connection coefficients and its derivatives one needs the spacelike geodesic deviation equation 



du 2 

It is convenient to write it as 



it 



cPnt 1 dT^ dn v ■ ■ 



Using deviation vectors 77^ = u8^ and t]q = 5q, we write ([39]) as follows: 



dV- 

-^Cu + 21 ;;,r + r* i* r ?u = //",„.cr'». 



"dP 0» A3' 



From Eq.(18) one obtains r A M „= <5°^ A ^- These conditions together with (f70|), ( [n] ) yield 











CP T^ni ( 



du 2 





-j2 



Finally, applying (JT^ - 17^) and ( p8[ ) we compute 



-\ Roao n»,\'.x>.x i - + i E% ,fc A".\ *.\ ' + ■ • • , 



4 rf - ^ + 1 Jh 3P + ± R° ijp n" k x i x^x k 

1 
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